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Abstract 

Taking into account spins, we classify all two-body non-relativistic Dark Matter 
annihilation channels to the allowed polarization states of Standard Model particles, 
computing the energy spectra of the stable final-state particles relevant for indirect 
DM detection. We study the DM masses, annihilation channels and cross sections 
that can reproduce the PAMELA indications of an e'^ excess consistently with the 
PAMELA p data and the ATIC/PPB-BETS e^ + e' data. From the PAMELA data 
alone, two solutions emerge: {i) either the DM particles that annihilate into W, Z, h 
must be heavier than about 10 TeV or (ii) the DM must annihilate only into leptons. 
Thus in both cases a DM particle compatible with the PAMELA excess seems to 
have quite unexpected properties. The solution {ii) implies a peak in the e"*" -|- e~ 
energy spectrum, which, indeed, seems to appear in the ATIC/PPB-BETS data 
around 700 GeV. If upcoming data from ATIC-4 and GLAST confirm this feature, 
this would point to a 0(1) TeV DM annihilating only into leptons. Otherwise the 
solution (i) would be favored. We comment on the implications of these results for 
DM models, direct DM detection and colliders as well as on the possibility of an 
astrophysical origin of the excess. 
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1 Introduction and highlights 



Cosmological observations imply that about 80% of the mass of the Universe is some 
unknown form of cold Dark Matter (DM) [1]. Presently the origin and nature of the DM 
particles, their mass, spin, couplings and other properties remain completely unknown. 
Among many possible cold DM candidates the most popular ones are the stable weakly 
interacting massive particles (WIMPs) which occur in many extensions of the Standard 
Model (SM), most notably in supersymmetry. If the DM WIMPs are thermal relics, their 
annihilation cross section must be crw ~ 3 ■ 10^^® cm^/sec, which indeed is typical of a 
weakly interacting TeV-scale particle. To discover the DM particles, experiments have 
searched for direct production of WIMPs in collider experiments, their scattering off the 
nuclei in terrestrial detectors as well as indirect signals of particles generated by the DM 
annihilation in the galactic halo [2]. 

The recently reported results by the PAMELA experiment [3] have the opportunity, 
if confirmed, of establishing a breakthrough in the cosmic antimatter searches. The 
PAMELA data show (i) a steep increase in the energy spectrum of the positron frac- 
tion, e"*"/ (e+ + e"), in cosmic rays above 10 GeV [3], compatibly with previous less certain 
hints from HEAT [4] and AMS-01 [5]; {ii) no excess in the p/p energy spectrum [6] com- 
pared with the predicted background; (Hi) at low energy, Ee+ < 10 GeV, the positron flux 
is presently suppressed by the solar magnetic polarity state A~ [7]. The positron excess 
observed by PAMELA is so big that it can produce observable spectral features in the 
total e"*" + e~ flux. Indeed, most recently the PPB-BETS balloon experiment reported an 
excess in the e"*" + e~ energy spectrum between 500-800 GeV [8] , confirming the similar 
earlier claim by the ATIC-2 balloon experiment [9]. These data-sets have not received 
much attention yet as probes of DM indirect detection, but, as we will see, they have 
the potential of being relevant to the issue. Of course one should be careful: the Monte 
Carlo simulations that such experiments need to tag and infer their energy have been 
tested only up to LEP energies; the excess is based on just a few data-points that are 
not cleanly consistent between ATIC-2 and the smaller PPB-BETS; emulsion chambers 
(EC) balloon experiments [10] do not show evidence for an excess, although they have 
larger uncertainties. The excess in both experiments seems to shows a cut-off at ener- 
gies just below 1 TeV. Those results will soon be tested by forthcoming ATIC-4 [11] and 
GLAST/Fermi data [12], that have bigger calorimeters able of containing a electromag- 
netic shower. All these excesses might be of course due to some astrophysical source, such 
as a nearby pulsar [13]. However, they also might provide the first clear evidence for the 
DM annihilations into SM particles. 

In this work we explore whether galactic DM annihilations can account for the observed 
excesses in the data. We perform a model independent analysis of the PAMELA positron 
and antiproton data with and without the electron plus positron data obtained from the 
balloon experiments. Our aim is to study if and how one can get model independent 
information on the DM mass, spin, interactions to the SM particles and on the DM 
annihilation cross section from the experimental data. Taking into account restrictions 
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that must be obeyed by non-relativistic DM annihilations, we identify the parameters 
that can be probed by the indirect signals. We study the non-relativistic DM annihilation 
cross sections into the set of all possible SM particle final states, DM DM SM SM, 
where 

SM = {e, /ii, fXR, tl,tii, Wl, Wt, Zl, Zt, h, q, b, t}, 

taking into account the allowed polarizations (Transverse, Longitudinal, Left, -Right). 
Using Monte Carlo tools, partly written by us, we try to keep the polarizations/helicities 
of the DM annihilation products and compute the energy spectra of the final state 
e^,p^, 7, *Fe,/i,r coming from their decays. Indirect signals of any DM model can be 
obtained by combining our decay spectra according to the branching ratios predicted by 
the model. The experimental results can be analyzed in terms of these phenomenological 
parameters, without committing to some specific theoretical scenario. 

Taking into account different models for the DM distribution and for the electron and 
proton propagation in the Galaxy, we perform fits of the predicted e"*"/ (e"*" + e~), p/p and 
(e"*" + e~) energy spectra to the PAMELA and balloon data. We allow for the possibility 
of enhanced annihilation cross sections and/or for the high density substructures in the 
DM halo, assuming an energy-independent boost factor B. For simplicity, we assume a 
common boost factor for the and for the and we take it as energy- independent. 
These are widely used assumptions, but they are not the only possibility [14]. One can 
easily infer what happens if these assumptions are relaxed. We thus study which DM 
masses and annihilation channels can best explain the experimental measurements. 

To illustrate our results, as well as to present the PAMELA and balloon e~^/ (e"^ + e~), 
(e++e~) and p/p data, we show in Fig. 1 three examples for the DM masses M = 150 GeV, 
1 TeV and 10 TeV that annihilate into W^W~ and //"""/i", as indicated in the figures. The 
first column presents the predictions for the positron fraction, the second for the total 
electron and positron flux and the third for the antiproton over proton flux. 

Many DM annihilation channels into SM primary channels can give a reasonably good 
fit to the PAMELA e"'"/(e+ + e~) excess alone, at least for some value of the DM mass 
picked in a wide range 60 GeV — > tens of TeV (first column of Fig. 1; see Fig. 8 for the 
full analysis). More precisely, DM annihilations into SM leptons can well fit for almost 
any mass in this range, while annihilations into quarks and Higgs bosons are disfavored at 
low masses because of the soft positron spectrum. Annihilations into gauge bosons work 
for light and heavy masses. 

The inclusion of the p data changes the issue. Because the DM annihilations to 
gauge bosons, Higgs bosons or quarks {W, Z, /i, g, 6, t) produce antiprotons, under our 
assumptions the PAMELA p/p data excludes light DM with sizable annihilation fraction 
into these channels (Fig. 1, the first row). However, for M > 10 TeV the antiproton 
spectrum becomes again consistent with present data, and the W^W~ seem to give the 
relatively best fit (Fig. 1, the third row). 

We find therefore that the PAMELA data single out two sets of solutions which can 
satisfy the observations compatibly with our standard astrophysical assumptions: (i) 
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DM with M = 150 GeV that annihilates into W^W 




Figure 1: Three examples of fits of (left), + e~ (center), p (right) data, for M = 
150 GeV (upper row, excluded by p), M = ITeV (middle row, favored by data), M = 
lOTeV (lower row, disfavored by the current e"^ + e~ excess). Galactic DM profiles and 
propagation models are varied to provide the best fit. See Sec. 4 for the discussion on the 
treatment of the uncertain astrophysical background. 



4 



a heavy ( > 10 TeV) DM that annihilates dominantly into W^W ; (ii) a DM that 
annihilates into SM leptons, with no strong preference for any DM mass. 

The inclusion of the balloon data qualitatively changes the sensitivity to the DM 
mass because the PAMELA excess predicts observable spectral features in the total flux 
(Fig. 1, the second column). It is striking that the same signal of the DM annihilation 
to leptonic channels which can explain the PAMELA data is also able to produce the 
apparent (e"*" + e~) excess in the PPB-BETS and ATIC data. Because the balloon data 
shows a sharp cut-off in the excess just below 1 TeV, the DM mass should be close to 

1 TeV, and all other but leptonic DM annihilation channels would be strongly disfavored 
or excluded. On the other hand, if the e"*" + e~ excess will not be confirmed, the e"*" + e~ 
data will strongly constrain or exclude the DM annihilations to leptons. Should this be 
the result of the upcoming ATIC-4 (or GLAST) experiment, the astrophysically favorite 
candidate for DM becomes a heavy ~10 TeV WIMP annihilating into SM gauge bosons. 

We comment on the implications of PAMELA and balloon data for supersymmetric 
DM candidates [15], as well as for some theoretically more exotic alternatives [16, 17]. We 
explore in a qualitative way the model-independent expectations for direct DM searches 
and for the production and detection of the DM particles at LHC and future lepton 
colliders. A few authors have discussed implications of the PAMELA e+/(e^ + e^) re- 
sults [18, 19, 20, 21]. While the analysis of the e+/(e+ + e") data presented in [20, 21] 
is consistent with our more systematic treatment, our final conclusions differ because we 
also include the important p and + e~ data. 

The paper is organized as follows. Section 2 discusses cosmological constraints on the 
DM annihilation cross section. In Section 3 we identify the DM annihilation channels and 
compute the spectra of SM final states. In Section 4 we discuss the particle propagation 
in the Galaxy. In Sections 5, 6 and 7 we perform fits to the PAMELA positron data, 
PAMELA positron plus antiproton data, and to the PAMELA and balloon data, 
respectively. We discuss the implications of our results for DM direct detection in Section 
8 and for collider phenomenology in Section 9. We conclude in Section 10. 

2 Cosmology and the total DM DM cross section 

The indirect signal rate is proportional to the total DM annihilation cross section. As well 
known, assuming that DM is a thermal relic, the measured DM cosmological abundance 
fixes the total DM annihilation rate [22]. Assuming that we can neglect co-annihilations 
with other new particles slightly heavier than DM, cosmology suggests 

(Tf ~ 3 X 10~^^cm^/sec around freeze-out, i.e. v ~ 0.2. (1) 

Astrophysical indirect DM signals are proportional to af at f ~ 10^^, which is the typical 
DM velocity (with respect to c) in our Galaxy, as determined by the escape velocity 
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Figure 2: Values of the DM annihilation cross sections suggested by the DM abundance. 

independently of the DM mass M. Fig. 2 illustrates the possible extrapolations of av 
from V ~ 0.2 down to f ~ IQ-^: 

a) av remains constant if, in the non-relativistic limit f <C 1, the cross section is domi- 
nated by s-wave annihilations. This leads to an indirect signal rate not much below 
the sensitivity of present experiments, that imply av < 10^^^ cm^/sec, depending on 
the branching ratios and on astrophysical issues. 

b) Alternatively, p-wave annihilations might dominate if s-wave annihilations are sup- 
pressed — in view of av oc f ^ this leads to a negligibly small rate of indirect DM 
signals. Case b) will therefore be ignored in our analysis. 

c) The most favorable possibility is that s-wave annihilations dominate and are en- 
hanced, in the non-relativistic limit, by the presence of Sommerfeld corrections, 
present in the SM if DM interacts with the W, Z vectors with a coupling g ^ g2 and 
is heavier than M>4'n'Mw/g'^ ~ 2TeV [23, 24]. In such a case av roughly grows 
as fmax/'W in the range Wmin < v < fmax where fmax ~ ^7^/4 and fmm ~ Mw/M [24]^ 

smaller Wmin is possible for specific values of M that lead to a zero-energy DM DM bound state; 
we estimate that DM annihilations in first structures [25] put the bound fmin ?^ 10^'^, compatibly with 
Minimal Dark Matter predictions, as shown in fig. 1 of [24]. 



6 




Figure 3: Contour plots of the Sommerfeld enhancement for one abelian massive vector 
with mass My and coupling a to DM (from [24]) and the resulting Sommerfeld reduction 
in the cosmological thermal DM abundance (we here assumed a DM mass M = 1 TeV). 

However details are model-dependent. The red line in Fig. 2 illustrates the Som- 
merfeld enhancement obtained for g = astrophysical signals are enhanced by a 
few orders of magnitude with respect to case a ) . 

We here proceed focussing on s-wave annihilations, and we assume the standard case a) 
to normalize our plots in Figs. 5, 6, 7. 

However, as we will see, the PAMELA and ATIC/PPB-BETS anomalies suggest that 
av and/or the electron boos factor is much larger than what is suggested by cosmol- 
ogy and/or astrophysics [14]. Co-annihilations only allow to gain (or loose) 0{1) factors. 
Possible resonances (e.g. a particle with mass very close to 2M ) may have a bigger effect 
but are strongly model dependent. If future data will confirm that such a sizable enhance- 
ment is needed, and disregarding possible resonances, one can pursue two possibilities: 
i) DM is not a thermal relic, so that it can have a av larger than what suggested by 
cosmology; ii) DM is a thermal relic, and av is Sommerfeld enhanced at low v. Within 
the SM this requires a DM mass heavier than a few TeV; in general one can assume that 
DM interacts with some extra light vector or scalar with mass My and coupling a to DM. 
Fig. 3 illustrates that a sizable Sommerfeld correction at low velocity v gives a relatively 
minor correction to the cosmological DM abundance. 
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3 Possible DM annihilations 



We study how Lorentz and gauge-invariance restrict the possible non-relativistic DM an- 
nihilations. We assume that the primary annihilation products contain two SM particles, 
so that the possible cases are 

W^W-, ZZ, Zh, hh, e+e", /^V", bb, ti, qq, (2) 

where q denotes any light quark, u, d, s, c. In view of the presumed neutrality of the 
DM particle we have not included final states containing photons. Bounds from direct 
detection experiments [26, 27, 28] suggest that the DM coupling to the Z is much smaller 
than g2, so that also the ZZ primary channel should perhaps be excluded from the list 
in eq. (2). 

Next, we need to consider the allowed polarizations of the SM primary particles in 
eq. (2), because the energy spectra of final-state e^,p^,7, ''i^i,fj.,T produced by their decays 
depend on their polarization. Since we assume that DM is an elementary weakly inter- 
acting particle, its spin is restricted to be or 1/2 or 1. In view of the spin-composition 
rules 

1®1 = 1, 2®2 = 1a©35, 3(E)3 = 1s®3a®5s, (3) 

the spin of the two-body non-relativistic s-wave DM DM state can be 0, 1 or 2. In the 
Majorana fermion case, the wave function must be asymmetric under exchange of the two 
DM particles, so that only the 1^ two-body state is allowed. In summary, for all DM 
spins the two-body DM DM state can have spin 0; higher spin is possible only in some 
cases. 

In order to be general and to avoid the lengthy expressions obtained from explicit 
evaluation of Feynman diagrams even in the simple DM models, we use tools from effective 
relativistic quantum field theory. Indeed, the non-relativistic nature of the DM DM 
two-body state is a relativistically invariant concept, and its annihilation products are 
quantistic and relativistic, so that these concepts are needed. In practice, all this amounts 
to say that the two body DM DM state is represented, depending on its spin, by a scalar 
field or a vector or ^^j^ if DM is a vector. This latter case is somewhat special, 
as it needs one extra gauge invariance, and is realized within Little Higgs models with 
T-parity as well as in models with extra dimensions.^ 

^ As the possible underlying motivations are not important for our purposes, we present a simplified 
version of its Lagrangian. In the simplest case the gauge group is abelian, e.g. a KK excitation of 
hypercharge. Taking into account this abelian factor of the SM, we introduce a gauge group G ~ 
U(l)i (8i 11(1)2 and three Wcyl fcrmions tp, ipi and -02 with charges (1,1), (—2,0), (0,-2) and a scalar 
H with U(l) charges (1, —1). Then the Lagrangian is given by the gauge-covariant kinetic terms for all 
fields, plus the {XitpiH + \2'>p2H*)ip Yukawa interactions. The vacuum expectation value of H breaks 
G to the SM U(l)y, leaving some particle massless (to be identified with the SM fermions and vectors, 
that acquire a mass when the SM gauge group is broken) and some other massive. Imposing a 1 <-> 
2 Z2 permutations symmetry enforces Ai = A2 such that the TJeavy and Light mass eigenstates are 
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In the light of these considerations, we can now hst the possible polarization states that 
give rise to observationally nonequivalent spectra of final SM particles, taking into account 
that DM annihilations cannot be observed event-by-event, but only on a statistical basis, 
after propagation in the Galaxy. 

1. DM annihilations into two massive vectors, W^W" or ZZ. We focus on the DM 
DM singlet state that exists for all possible DM spins. It can couple to SM 
vectors in 3 possible ways, described by the following effective Lagrangians: 

la) The effective interactions ^F^^e^^p„Fp„ and ^-F^^ both give rise to annihilation 
into vectors with transverse polarizations (helicity correlations are different in 
the two cases, but cannot be observed). Transversely polarized massive vectors 
decay in two light fermions with energy E = x M as: 

dN/dcose = 3{l + cos^e)/8 i.e. dN / dx = 3{l - 2x + 3x^) /2, (4) 

where distributions are normalized to = 1, ^ is the angle between the mo- 
menta of the vector and of the secondary fermion and the last expression holds 
for M > Mw 

lb) The possible effective interaction gives vectors with longitudinal polar- 
izations. Within an SU(2)l invariant theory, such effective couplings arises 
when DM couples to the components of the Higgs doublet that get eaten by 
the W, Z vectors. A longitudinal vector decays as 

dN / d cos 9 = 3{l-cos'^ 9)/ 4 i.e. dN/dx = 6x{l - x). (5) 

Fig. 4a shows the energy spectra of secondary fermions, in the limit M ^ Mw- 

2. DM annihilations into the Higgs bosons. We can again focus on ^, so that the 
effective interaction ^h'^ gives DM annihilations into two Higgs bosons. Since they 
have no spin, there are no polarization issues. Precision Electro Weak data suggest 
a light Higgs boson mass: we here assume that the Higgs mass is 115 GeV so that 
the Higgs dominantly decays into b and r. (Higgs decay modes into heavy vectors 
would be instead dominant if mh > 2Mw)- DM annihilations into Zh will not be 
considered, as they are are essentially given by the average of the ZlZl and hh 
channels. 

Wh = Wi — W2, Wl = Wi + W2, i'L — V"! — "02, i^H — V'l + "^2 and the final Lagrangian has the 
form C ^ H^Wfj + Hijjipn + WhiPh^l + WLiip"^ + V'i + V'lf)- The Hghtest among the heavy vector 
Wh or fermion ^pn is stable thanks to the Z2 parity. A more complicated construction can be applied 
to the non-abelian factors of the SM allowing to get light chiral fermions, and again one finds potential 
vector or fermion DM candidates. UV completions the little-Higgs models usually render the permutation 
symmetry anomalous [29], ruining the stability of DM. In the Universal Extra Dimension context, the 
problem is that truly universal extra dimensions do not give chiral fermions, so that one needs to add 
extra structures such as boundaries: divergent quantum corrections generate operators there localized, 
ruining the 'universality'. 
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Figure 4: Energy spectra of fermions f produced by decays of a heavy vector (e.g. f = e 
from W e~^i>e) or of a heavy fermion (e.g. f = e from /i+ — > e^VeV^) from polarized 
vector and fermion decays, in the limit M ^ Mw,z ^ ruf. The dashed lines show the 
result for the unpolarized case. 

3. DM annihilations to SM fermions, DM DM — > //, with mass rrif. 

3a) The DM two-body singlet state 3^ can only couple as ^/i/j^ + h.c. = ^^^/\E'/, 
where Dirac and Weyl notations for fermions are related by \E'/ = (/l, fn). If 
mj is not negligible, the operator can have a complex coefficient with physical 
meaning; the imaginary part gives rise to the 3^ f'^^^ operator. Anyhow, the 
h.c. implies that on average the produced fermions have zero helicity. Further- 
more, the Weyl notation makes clear that, unless DM is coupled to both fi 
and /i?, the coefficient of this operator is suppressed by rrif/M. 

3b) We next need to consider the vector two-body DM state, present if DM 
is a vector or a Dirac fermion. It can couple to fermions in two different 
ways: ^^[/l7m/l] = n^^PL"^] or as %[fR^M = ^I^^^Pr^ The 
first (second) operator means that / fermions have negative (positive) helicity. 
To see how helicity affects the energy distribution of final-state fermions, let 
us consider weak decays / — >• f'f'f" such as /z"*" — >• p^e'^Pe in the limit of 
mf/,mf",mfi'/ -C m/ -C Mw- In view of Fierz identities, /' and /" have the 
same energy distribution: 

dN/dx\L = 2{1 - xf {I + 2x), dN/dx\R = 4:{l-x^)/3. (6) 
We can consider the two limiting cases of L and R spectra because, as plotted in 
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Fig. 4b, their average is the unpolarized spectrum of point 3a). 



Having precisely identified the primary two-body states that can be produced in DM 
annihilations, we next need to study how the spectra of final-state e,p, v can be computed 
using MonteCarlo codes such as Pythia [30], that mostly neglect polarizations. We 
computed polarized W, Z decays using our own MC routine, and decayed the resulting 
fermions using PythiaS. For the /i spectra, we used analytic results. For the r spectra 
we used Tauola [31], interfaced with Pythia6. It would be interesting to repeat such 
computations with other MC tools that employ different hadronization models, in order 
to get an estimate of the theoretical uncertainty. 

In the simplest case of DM annihilations into e+e^, we have taken into account radi- 
ation so that electrons and positrons have a smooth peak at a; < 1 rather than a line at 
a; = 1, and the photon spectrum is described by the Altarelli-Parisi splitting function 



dN^ a 
dx 2tt 



2 



l + ln( — ^ 1-x 

mi 



1 + (1 -x) 



X 



(7) 



These are the most interesting annihilation channels that will emerge from our later anal- 
ysis, and, especially in the n case, polarizations are relevant. For quarks, the polarization 
issue should be irrelevant compared to hadronization. A possible exception is the t quark. 

In this way we obtained the e and p spectra at DM annihilations, plotted in Fig. 5 
and Fig. 6. These spectra will undergo deformations due to the propagation effects of the 
charged particles in the galactic halo, to be discussed next. 



4 Propagation in the Galaxy and the astrophysical 
background 

With the computed primary spectra in hand, we next need to consider DM annihilations 
in the Galaxy, and propagation of the resulting and p up to the solar system. This is 
done proceeding as in [32, 33, 34], so we just briefly summarize the main points. Three 
different halo models are considered: isothermal, NFW and Moore. For each one of them, 
three different sets of parameters for propagation (different for and p) are considered: 
MIN, MED and MAX. The main physical parameter that discriminates between these 
conflgurations is the thickness of the cylinder, centered on the galactic plane, inside which 
charged particles such as and p diffuse. The remaining parameters, such as the diffusion 
coefficients for e^,p, the energy losses for and the galactic convective wind speed for 
p, are fixed in order to reproduce the measured spectra of Cosmic Ray nuclei. 

In line of principle, the astrophysical background spectra should be computed together 
with Cosmic Ray nuclei, such that uncertainties are properly correlated. We instead as- 
sume that the e~^,e~ ,p background spectra can be freely renormalized, and have indepen- 
dent ±0.05 uncertainties in their energy slope (namely, the central values used in [34] are 
multiplied by AE^ and the resulting is minimized with respect to p = ± 0.05 and 
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Figure 5: Energy spectra of the positron fraction e+/(e^ + e^), as produced by DM an- 
nihilations into the various possible channels. The DM masses are {50,100,200,500,- 
1000, 2000, 5000, 10000, 20000} GeV and can be inferred from the energies at which the 
spectra drop to zero. In this figure we take a NFW halo, the MED propagation model, 
boost factor = 1 and the cross section av = 3 ■ 10^^^ cm'^/sec naively suggested by 
cosmology. For illustration, the upper plots show the expected astrophysical background 
(upper line) and the data, that can be fitted increasing Be ■ av. The dashed curves in the 
/i, r, W, Z plots represent fifi, tr, Wl, Zl polarizations, while the continuous lines represent 
fJ-L, tl, Wt, Zt polarizations. 
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Figure 6: The same as in Fig. 5 for the energy spectrum of the p/p ratio. 



to A). This mimics the main uncertainties in astrophysical backgrounds'^, produced by 
the Fermi mechanism of acceleration, that typically generates power-law spectra (up to 
some cut-off) but does not predict its coefficient. The positron fraction below 20 GeV has 
been measured very precisely, but is strongly affected by solar modulation, which is not 
understood well enough to allow a useful theoretical computation. Therefore we assume 
that each data-point for the positron fraction has an independent uncertainty given by 
the spread between different solar modulations, as plotted in [7]. This amounts to a ±6% 
uncertainty at 10 GeV and ±30% at 1 GeV, and in practice means that low energy data 
do not provide useful constraints on Dark Matter annihilations. 

In performing the fit, we smoothly scan over the intermediate propagation configura- 
tions and halo models, within the boundaries described above. The MED configuration 

^We checked that this procedure reproduces reasonably well the uncertainty bands reported by more 
detailed analysis, see e.g. Fig. 7 of [35] for the case of antiprotons and Fig. 4 of [36] for positrons. It seems 
to us unlikely that uncertainties in the electron astrophysical background (PAMELA is measuring it more 
accurately than ATIC [9]) can produce the rising positron fraction apparent in the PAMELA data. 
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Figure 7: Comparison of the electron (left) and proton (center) fractions and photon 
(right) fluxes produced by possible DM annihilation channels, for M = 1 TeV. 



is sometimes considered as favored, but we do not attach a statistical meaning to this 
sentence. 

Marginahzations over nuisance parameters and other statistical operations are per- 
formed as described in Appendix B of [37]. We will show plots of the &s a function of 
the DM mass: an interval at n standard deviations corresponds (in Gaussian approxima- 
tion) to < Xmin + "-^5 irrespectively of the number of data points. We will not report 
the value of x^/dof as it is a poor statistical indicator; furthermore the number of dof 
is not a well-defined quantity when (as in the present case) data-points with accuracies 
much smaller than astrophysical uncertainties are effectively irrelevant. 

5 PAMELA positron data 

We start our data analysis considering only the PAMELA e+/(e"'" + e~) observations (16 
data points) [3]. 

Taking into account the DM distribution and positron propagation effects in the 
Galaxy, the energy spectra of the positron fraction originating from different DM an- 
nihilation channels is plotted in the left panel of Fig. 7 for the DM mass M = 1 TeV. 
As expected, the most energetic positrons come from the pure leptonic channels and the 
softest spectra are produced in quark annihilation channels. 

Fitting data as described in the previous section. Fig. 8 shows how well the possible 
DM annihilations into two SM particles can fit the PAMELA positron excess. Fig. 9 
shows the boost factor (with respect to the cross section suggested by cosmology, 
av = 3 10"^^ cmVsec) and ■ av that best fits the PAMELA excess. We see that DM 
annihilations into e, fi, r, W can reasonably well reproduce the data for any DM mass. 
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Figure 8: A fit of different DM annihilation channels to the PAMELA positron fraction 
data only. The labels on each curve indicate the primary annihilation channel. 



while annihilations into Z, t, q, b, h give a good fit for DM heavier than about 1 TeV. It 
is perhaps interesting to note that, contrary to what commonly thought, the spectrum 
from W^W~ annihilations is not too flat to give a good fit of the quite steep PAMELA 
rise. At small masses (see e.g. the upper-left panel of Fig. 1) a MIN configuration of the 
propagation parameters (and a proper variation of the background curve within the limits 
considered above) allows to fit the data. At large DM masses (see e.g. the lower-left panel 
of Fig. 1) the low-x portion of the primary spectrum is steep enough to do the job (as 
usual, X = E/Mdm)- 

It is interesting to study if the PAMELA excess can be produced by supersymmet- 
ric DM. Already at this stage one can see that it is necessary to give up some usual 
assumptions, either the naturalness of supersymmetry (invoking multi-TeV DM) and/or 
DM astrophysics (invoking very large boost factors, as in [18]) and/or that the DM is a 
thermal relic (invoking some non-thermal DM cosmological production mechanism). In 
this last case one interesting possibility emerges [38]. A pure Wino triplet annihilates 
predominantly into transverse W^W~ with 

27rM2(2-M2,/M2)2' 

that for M ^ 100 GeV is precisely the value suggested by the PAMELA e"*" excess, that 
would therefore be naturally produced with the boost factor of order unity. 

We now consider the PAMELA p data, showing that this latter and many other 
possibilities are excluded because they lead to large unseen p excess, as already illustrated 
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DM mass in GeV 

Figure 9: Values of Bg-av (right axis) and of the boost factor (left axis, for the standard 
av = 3 10~^^ cm^/secj suggested by the PAMELA excess. 

in Fig. 1. 

6 PAMELA positron and anti-proton data 

We address now the implications of the PAMELA p/p data [6] (17 data points). The 
middle panel of Fig. 7 shows the p/p fractions produced by the various primary channels. 
One can see that there are two basic features: i) annihilations into leptons do not give rise 
to protons; ii) all other channels give rather similar proton energy spectra. Comparison 
with the positron spectra reveals that the proton energy spectrum is softer than the 
positron one produced in the leptonic and gauge boson annihilation channels. 
The results of the fits are shown^ in Fig. 10. 

Consistently with what anticipated, since no excess seems present in the p/p ratio, 
annihilation into leptons are not constrained as they do not produce antiprotons. On 
the contrary, all other annihilations into quarks, vector and Higgs bosons are significantly 
constrained, and allowed only if the DM particle is heavier than almost 10 TeV. Only in 
such a case the proton excess lays at energies above those explored currently by PAMELA, 
while the low energy proton spectrum is consistent with the background (see Fig. 1 for 
illustration). The bound dominantly comes from high energy data points where the solar 

^In order to study how the fit of PAMELA data changes adding p/p data we assume equal boost factors 
and propagation models for positrons and antiprotons, as discussed in Sec. 4. The second assumption 
should be a good approximation, while the first one can easily fail when boost factors are very large, 
allowing to relax all our results. 
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Figure 10: Combined fit of different DM annihilation channels to the PAMELA positron 
and PAMELA anti-proton data, assuming equal boost factors and propagation model. 



modulation is negligible. 

The implications of the complementarity of PAMELA e^/{e~^ + e~) and p/p data on 
constraining new physics are evident. The light Wino considered in the previous section 
is excluded as a DM candidate because its annihilations would induce a large unobserved 
antiproton excess. Let us consider models where DM is a vector, the lightest Kaluza-Klein 
(KK) state or some other heavy replica of the photon or of the hypercharge boson or of 
the Z. They do annihilate directly into leptons, but the direct annihilation rates into 
quarks are of the same magnitude or larger, as quarks have 3 colors and similar charges. 
For example the heavy KK of the hypercharge annihilates into fermions // with s-wave 
cross section 

Such a p excess cannot escape detection at PAMELA unless the DM particles have multi- 
TeV masses. The same conclusion applies to the extra dimensional models in which the 
DM is the lightest KK state of neutrino as those annihilate predominantly to WW final 
states. Similarly the DM candidates in the Little Higgs models must be very heavy as they 
always have annihilation channels producing antiproton excess. Therefore, the PAMELA 
data alone put very strong constraints on the mass of DM that have significant fraction 
of annihilations into quark or gauge boson channels. 
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DM mass in GeV 

Figure 11: Combined fit of PAMELA positron and balloon + e~ data. 

7 PAMELA positron and balloon e+ + e~ data 

Various balloon experiments measured the spectrum of + cosmic rays at energies 
higher than those where PAMELA can measure separately the e~ and spectra. If the 
e"*"/ (e"*" + e~) excess present in PAMELA data continues to rise in the same way at higher 
energies, the positron fraction becomes of order unity at a few hundred GeV, giving rise 
to an excess in the e~ + e~ spectrum. The most recent and precise measurements by 
ATIC-2 [9] and PPB-BETS [8] indeed show a hint of an excess at such energies, which is 
not visible in the less precise EC data [10]. The experimental situation is therefore quite 
open at these high energies, due also to the intrinsic difficulty of the measurements and 
the analysis. The ATIC-4 [11] experiment already collected new data, and GLAST [12] is 
collecting new precise data, so that the possible structure in the e"*" + e~ spectrum should 
soon be clarified. Keeping this in mind, we explore the impact of the currently available 
results. We fit ATIC-2, PPB-BETS and EC data, for a total of 37 data points. 

As the DM annihilations necessarily should affect also the total + e~ energy spec- 
trum, various issues deserve to be studied. 1) Are DM annihilations compatible with this 
hint? 2) Are DM annihilations compatible with this hint and with the PAMELA excess? 
3) Even if this e~ + e~ hint will not be confirmed, do e~ + e~ data already constrain DM 
annihilations? The answer to all these questions is yes. 

A comparison of the energy spectra produced by the various annihilation channels 
shown in Fig. 7a indicates that only DM annihilations into e, /z or r can reproduce the 
e"*" + e~ excess for a DM mass around 1 TeV. In this case one obtains an excess in the 
positron fraction quantitatively similar to the one present in the PAMELA data. In 
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Fig. 11 we present the combined fit to the lepton data. 

• DM annihilations into n seem to give the optimal energy spectrum and the best fit, 
as illustrated in the middle row of Fig. 1. Annihilations into e (r) give a slightly 
poorer fit, because of a too (not enough) steep spectrum. The DM mass has to 
be around 1 TeV. Indeed, a heavier multi-TeV DM that annihilates into e, /i, r is 
excluded, as a solution to the PAMELA excess, by + data. A lighter DM that 
fits the PAMELA excess instead gives a > 20% drop in the + e~ spectrum, not 
present in the ATIC-2 data: this is illustrated by the upper row of Fig. 1. 

• All other non-leptonic channels give an + e~ excess which is too smooth to 
reproduce the hint in e"*" + e~ data, as illustrated in the bottom row of Fig. 1: 
should the excess in the + e~ spectrum not be confirmed, multi-TeV DM with 
dominant non-leptonic annihilations will be able to reasonably agree with a smooth 
e"*" + e~ spectrum, while still producing the PAMELA e^/{e~^ + e~) excess. 

So, as the combined fit to all existing data prefers 1 TeV DM which annihilates only 
to leptons, all the popular DM candidates discussed so far in this paper are disfavored. 
If ATIC-4 will confirm the + e~ anomaly, DM particles must have quite unexpected 
properties. 



8 Implications for DM direct detection 

We explore the connection between the cross section Uann'^^ ~ 3 ■ 10^^^ cm^/sec ■ (M/ TeV)^ 
suggested by the PAMELA indirect signal and the cross section o"dir probed by direct 
detection experiments. Rather than studying precise but model-dependent implications, 
we try to delineate the qualitative generic pattern. Depending on the DM annihilation 
channel, there are three different possibilities. Taking into account a possibly large boost 
factor Be and/or Sommerfeld enhancement S, we estimate: 

i) DM annihilation into particles contained in nuclei, i.e. into quarks. The qualitative 
connection between the cross sections is 



(9) 



m%aannV 10 ^^cm^ 

AP BeS BeS 

although specific couplings can give one extra mj^/M'^ suppression. The exper- 
imental bound on the spin-independent direct cross section is at the 10~^^ cm^ 
level [26, 27]. In the light of large boost factors B^. and/or Sommerfeld correction 5* 
required by PAMELA, eq. (9) can be compatible with the present bounds. 

ii) DM annihilations into unstable SM particles not contained in normal matter, such 
as W, Z, h, fi, T. Then a direct cross section is generated at loop level, so that 
one expects a adh smaller than in eq. (9) by one or two loop factors ~ (47r)^^. 



19 



Furthermore, the mass scale at the denominator can be rather than M^. The 
connection is therefore model dependent: e.g. the Minimal Dark Matter quintuplet 
predicts (Jdir ~ 10~'^*^cm^ [39]. 

iii) DM annihilations into electrons. In order to observe the DM elastic scattering on 
electrons bounded in atoms, DM DM e~, one needs collisions with the most 

bound atomic electrons, that have quasi-relativistic momentum, p ~ rrie, such that 
the recoil energy AE ~ vp is above threshold [40] . The qualitative connection is 

mla^nnV 10"^^ cm^ 

The experimental apparatus must be able to detect electrons and must have low 
energy threshold for the electron detection. Such an experiment is, for example, 
DAMA [28]. In order to explain the DAMA anomaly, the DM e~ cross section 
should be ddir ~ 10"^^ cm2(M/ TeV). Thus, in this scenario, the PAMELA data 
points to very small direct cross section compared to the sensitivities of the present 
experiments. 

Therefore DM motivated by the PAMELA anomaly is easily compatible with all the 
constraints from direct detection experiments. 



9 Implications for collider experiments 

We finally briefly discuss the implications for collider searches of the DM possibilities that 
our analysis has individuated as favored by the current data. 

DM masses of (9(10) TeV are of course out of the reach of the LHC DM particles 
with ~1 TeV mass and weak interactions are at the limit of what can be detected by LHC 
experiments. DM particles with ~1 TeV mass which only couple to leptons would escape 
detection at LHC [41], unless accompanied by extra slightly heavier colored particles. In 
the latter scenario the experimental signature is the usual missing transverse momentum 
carried away by DM in the cascade decays of the new model-dependent colored particles. 
Such a DM cannot be produced in the ILC either, because of kinematical limitation 
on the center of mass energy y/s = 0.5 — 1 TeV. At future hypothetical e"'"e~ colliders 
like CERN CLIC or in the fi~^fi~ colliders with ^/s > 2 TeV, the direct production of 
DM particles is feasible. In particular, the PAMELA results on the annihilation cross 
section DM DM— imply a lower bound on the DM pair production cross section 
a{i^i~ — s> DM DM) in those colliders. Even if the required large boost factor B is of 
astrophysical or cosmological origin, the production cross section of order 1 pb is large 
enough to observe such a DM at lepton colliders. 
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10 Discussion and conclusions 



We studied if DM annihilations into two-body SM states can reproduce the features that 
seem present in the energy spectra of e"^ , e~ , p cosmic rays recently measured at energies 
between about 10 GeV and 3 TeV by the PAMELA [3], ATIC-2 [9] and PPB-BETS [8] 
experiments. 

We found that: 

• Considering the PAMELA positron data alone, the observed positron excess can be 
well fitted by DM annihilations into W, Z, e, fi, r with any DM mass, and by DM 
annihilations into q, b, t, h with a multi-TeV DM mass. 

• Adding the PAMELA anti-proton data suggests that the positron excess can be 
fitted by DM annihilations into W, Z, h (annihilations into q, b, t give poorer fits) 
only with a multi-TeV DM mass, unless boost factors or propagation strongly dif- 
ferentiate between e"*" and . DM annihilations into leptons are still viable for any 
DM mass. 

• Adding the balloon e"*" + e~ data similarly suggests that DM annihilations into 
W, Z, h, b, q, t can fit the PAMELA excess only with a multi-TeV DM mass. An- 
nihilations into e,fi,T predict a feature in the e"*" + e~ spectrum that, for a DM 
mass around one TeV, is compatible with the hint present in ATIC-2 and PPB- 
BETS data. A TeV mass seems too heavy for explaining the DAMA anomaly [28] 
compatibly with CDMS bounds in terms of inelastic DM scatterings [46]. 

• The annihilation cross section suggested by the PAMELA data is a few orders of 
magnitude larger than what naturally suggested by the cosmological abundance, 
unless DM formed sub-halos. This enhancement could be due to a Sommerfeld non- 
relativistic enhancement of the annihilation cross section. Alternatively, this may 
indicate for non-thermal production of DM. 

• Present data start to be sensitive to the polarization state of the DM annihilation 
products. 

In the light of our results two solutions emerge for the DM. 

1. The first possibility is a heavy DM, M> lOTeV, that annihilates into W^W~ or 
hh (an example of a DM candidate of this sort is provided by the model of [39]). 
This possibility is illustrated by the lower row in Fig. 1, and will be tested by future 
p data. 

2. The second, more exciting, possibility is that the PAMELA excess has also been 
seen in + e~ data, by ATIC and PPB-BETS. The forthcoming ATIC-4 and 
GLAST data should soon clarify this issue. In such a case the best fit is obtained 
for M ^ 1 TeV with DM annihilating into ^~ , and good fits are also obtained 



21 



for M ^ 800 GeV if DM annihilates into e+e", or M ^ 2TeV if DM annihilates 
into r+r^. This is illustrated by the middle row in Fig. 1. The needed 'boost times 
cross section' is SgCrt' ~ 3 10^^^ cm^/sec. 

In both DM particle compatible with the PAMELA anomaly has therefore unex- 

pected properties. In general, p and + 6" data seem to indicate that DM masses around 
100 GeV cannot fit the PAMELA anomaly, as illustrated by the fit in the upper row of 
Fig. 1. The favorite DM candidate, the supersymmetric neutralino [15], cannot annihilate 
directly into light leptons with a large cross section [18]. A multi-TeV Wino is a possible 
supersymmetric candidate, but supersymmetry would not naturally explain the breaking 
of the electroweak symmetry. The vector or Dirac fermion DM candidates appearing in 
the models of extra dimensions [42, 43, 44] and in the Little Higgs models [45] can anni- 
hilate directly into leptons. However, at the same time they also annihilate into quarks 
and/or into gauge and Higgs bosons. Therefore the masses of such DM candidates are 
also pushed to the multi-TeV region by the PAMELA antiproton data. Similarly, a LSP 
right-handed sneutrino with a sizable vrLH Yukawa coupling annihilates into both L and 
H (and consequently into longitudinal Z, that contain the Goldstone components of 
the Higgs doublet if), unless lighter than M\y. None of these theoretically motivated DM 
candidates mentioned here can explain the ATIC/PPB-BETS excess. 

A model of a thermal DM candidate that can explain both the PAMELA and the 
ATIC/PPB-BETS excesses with an order one boost factor can be constructed as follows. 
We assume that lepton fiavor — L,- is a spontaneously broken anomaly-free U(l)i;, 
gauge symmetry, with gauge coupling ai ~ 1/50 and vector mass ^ Mz- The DM 
candidate is a Dirac fermion with mass M ^ 1.5 TeV, charge ~ 2 under the extra 
U(l), and no SM gauge interactions. Its tree-level annihilation cross section relevant 
for cosmology, av = vr(3g|^ + qj^)a\/2M'^ = 3 10~^^cm^/sec at f ~ 0.2, is enhanced in 
astrophysics, up to av < 10^"^^ cm^/sec at f ~ 10^'^, by the Sommerfeld correction due to 
the extra U(l), as computed in Fig. la of [24]. The photon spectrum resulting from these 
DM annihilations into fi and r is hard, as can be seen in the right panel of Fig. 7. The 
correction to the magnetic moment of the muon, = +5f|m^/127r^M|, reproduces the 
observed discrepancy with respect to the SM, (5a^ = (3 ± 1)10~^ [47]. The new vector 
affects precision data as Re = Y{Z i^i-)lY{Z e+e ) = 1 + 3aL/(Mi/Mz)/47r 
with /(O) = 1; this is compatible with the LEP measurements, = 1.001 ± 0.003 and 
Rr = 1.002 ± 0.003 [48]. Furthermore one loop-effects generate a kinetic mixing of the 
new vector with the photon, 6 = e5'Lln(mT-/m^)/67r^ ~ 0.005, that gives rise to a cross 
section for direct detection, asi = 4:7iq1aLcym%9'^ /My ~ 10~'^^cm^, around the present 
bound. Reducing My by a factor of few leads to a too large asi, and reducing My below 
about 1 GeV also to unseen DM annihilations in first structures [25] and to coUisional 
DM with elastic a ~ g^M'^ /AtiM^. 

Finally, a third, less exciting possibility is that the features that appear in the e"*" and 
e"*" + e~ spectra are not due to DM annihilations but to some astrophysical object, such 
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as a pulsar. Indeed it is expected that pulsars produce a power law spectrum (with index 
p around 1.5-2) of mostly electron-positron pairs, with a cut-off. The black line in Fig. 11 
shows that a flux $e- = $,+ = A ■ E-^e"^/*^ with p ~ 2 and cut-off at M ~ 300 GeV 
does fit the e~^,e~ data, giving a peak less sharp than what hinted by ATIC-2 data. It 
is remarkable that pulsars naturally give an excess only in e^, with the observed slope 
p. While the total energy and total power emitted by a pulsar can be reliably computed 
assuming that they are respectively dominated by rotation and magnetic dipole radiation, 
the factors A and M cannot be computed without a detailed model of the acceleration 
mechanism. Furthermore, while the known nearby pulsars Geminga or B0656+14 are 
the main candidates, unknown pulsars can contribute. According to [13], young nearby 
pulsars tend to give a cut-off M at energies higher than what suggested by ATIC/PPB- 
BETS, while older pulsars tend to give a too low flux, as had time to diffuse away. 
An old intense pulsar no longer emits 7, so that it is difficult to distinguish it from DM 
annihilating into leptons; hopefully angular anisotropies of the e"*", e~ fluxes will help [13]. 
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